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. We apply a new general method of anholonomic frames with associated nonlinear connection 

structure to construct new classes of exact solutions of Einstein-Dirac equations in five dimensional 
(5D) gravity. Such solutions are parametrized by off-diagonal metrics in coordinate (holonomic) 
bases, or, equivalently, by diagonal metrics given with respect to some anholonomic frames (pentads, 
or funfbein, satisfing corresponding constraint relations). We consider two possibilities of general- 
ly-^ , ization of the Taub NUT metric in order to obtain vacuum solutions of 5D Einsitein equations with 
effective renormalization of constants (by higher dimension anholonomic gravitational interactions) 
having distinguished anisotropics on an angular parameter or on extra dimension coordinate. The 
. constructions are extended to solutions describing self-consistent propagations of 3D Dirac wave 
^ ' packets in 5D anisotropic Taub NUT spacetimes. We show that such anisotropic configurations of 
spinor matter can induce gravitational 3D solitons being solutions of Kadomtsev-Petviashvili or of 
sine-Gordon equations. 

cn 
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I. INTRODUCTION 

D . Recently one has proposed a new method of construction of exact solutions of the Einstein equations on (pseudo) 
■ Riemannian spaces of three, four and extra dimensions (in brief, 3D, 4D,...), by applying the formahsm of anholo- 
nomic moving frames There were constructed static solutions for black holes / tori, soliton-dilaton systems and 
wormhole / flux tube configurations and for anisotropic generalizations of the Taub NUT metric Q ; all such solutions 
being, in general, with generic local anisotropy. The method was elaborated following the geometry of anholonomic 
^ frame (super) bundles and associated nonlinear connections (in brief, N-connection) |^,^ which has a number of 
" " ' applications in generalized Finsler and Lagrange geometry, anholonomic spinor geometry, (super) gravity and strings 
with anisotropic (anholonomic) frame structures. 

In this paper we restrict our considerations for the 5D Einstein gravity. In this case the N-connection coefficients 
are defined by some particular parametrizations of funfbein, or pentadic, coefficients defining a frame structure on 
(pseudo) Riemannian spacetime and describing a gravitational and matter field dynamics with mixed holonomic 
(unconstrained) and anholonomic (constrained) variables. We emphasize that the Einstein gravity theory in arbitrary 
dimensions can be equivalently formulated with respect to both holonomic (coordinate) and anholonomic frames. In 
the anholonomic cases the rules of partial and covariant derivation are modified by some pentad transforms. The 
point is to find such values of the anholonomic frame (and associated N-connection) coefiicients when the metric is 
diagonalized and the Einstein equations are written in a simplified form admitting exact solutions. 

The class of new exact solutions of vacuum Einstein equations describing anisotropic Taub NUT like spacetimes 
1^ is defined by off-diagonal metrics if they are given with respect to usual coordinate bases. Such metrics can be 
anholonomically transformed into diagonal ones with coefficients being very similar to the coefficients of the isotropic 
Taub NUT solution but having additional dependencies on the 5th coordinate and angular parameters. 

We shall use the term locally anisotropic (spacetime) space (in brief, anisotropic space) for a (pseudo) Riemannian 
space provided with an anholonomic frame structure induced by a procedure of anholonomic diagonalization of a 
off-diagonal metric. 
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The Hawking's |^ suggestion that the EucHdean Taub-NUT metric might give rise to the gravitational analogue 
of the Yang-Mills instanton holds true on anisotropic spaces but in this case both the metric and instanton have 
some anisotropically renormalized parameters being of higher dimension gravitational vacuum polarization origin. 
The anisotropic Euclidean Taub-NUT metric also satisfies the vacuum Einstein's equations with zero cosmological 
constant when the spherical symmetry is deformed, for instance, into ellipsoidal or even toroidal configuration. Such 
anisotropic Taub-NUT metrics can be used for generation of deformations of the space part of the line element defining 
an anisotropic modification of the Kaluza-Klein monopole solutions proposed by Gross and Perry ||^ and Sorkin [0| . 

In the long-distance limit, neglecting radiation, the relative motion of two such anisotropic monopoles can be also 
described by geodesic motions, like in Ref. but these motions are some anholonomic ones with associated nonlinear 
connection structure and effective torsion induced by the anholonomy of the systems of reference used for modelling 
anisotropics. The torsion and N~connection corrections vanish if the geometrical objects are transferred with respect 
to holonomic (coordinate) frames. 

From the mathematical point of view, the new anholonomic geometry of anisotropic Taub-NUT spaces is also 
very interesting. In the locally isotropic Taub-NUT geometry there are four Killing- Yano tensors |9j. Three of 
them form a complex structure realizing the quaternionic algebra and the Taub-NUT manifold is hyper-Kahler. In 
addition to such three vector-like Killing- Yano tensors, there is a scalar one which exists by virtue of the metric 
being of class D, according to Petrov's classification. Anisotropic deformations of metrics to off-diagonal components 
introduce substantial changes in the geometrical picture. Nevertheless, working with respect to anholonomic frames 
with associated nonlinear connection structure the basic properties and relations, even being anisotropically modified, 
are preserved and transformed to similar ones for deformed symmetries Q. 

The Schrodinger quantum modes in the Euclidean Taub-NUT geometry were analyzed using algebraic and analytical 
methods I^Jl^ . The Dirac equation was studied in such locally isotropic curved backgrounds . One of the aims 
of this paper is to prove that this approach can be developed as to include into consideration anisotropic Taub-NUT 
backgrounds in the context of the standard relativistic gauge-invariant theory of the Dirac field. 

The purpose of the present work is to develop a general S'0(4, 1) gauge-invariant theory of the Dirac fermions |T^ ] 
which can be considered for locally anisotropic spaces, for instance, in the external field of the Kaluza-Klein monopole 
|pj| which is anisotropically deformed. 

Our goal is also to point out new features of the Einstein theory in higher dimension spacctimc when the locally 
anisotropic properties, induced by anholonomic constraints and extra dimension gravity, are emphasized. We shall 
analyze such effects by constructing new classes of exact solutions of the Einstein-Dirac equations defining 3D soliton- 
spinor configurations propagating self-consistently in an anisotropic 5D Taub NUT spacetime. 

We note that in this paper the 5D spacetime is modeled as a direct time extension of a 4D Riemannian space 
provided with a corresponding spinor structure, i. e. our spinor constructions are not defined by some Clifford algebra 
associated to a 5D bilinear form but, for simplicity, they are considered to be extended from a spinor geometry defined 
for a 4D Riemannian space. 

We start in Section II with an introduction in the Einstein-Dirac theory formulated with respect to anholonomic 
frames with associated nonlinear connection structure. We write down the Einstein and Dirac equations for some 
classes of metric ansatz which can be diagonalized via anholonomic transforms to a locally anisotropic basis. In 
Section III we construct exact solutions for 3D Dirac wave packets propagating in anisotropic backgrounds. In 
Section IV we outline two classes of 5D locally anisotropic solutions of vacuum Einstein equations generalizing the 
well known Taub NUT metric for the cases of anisotropic (angular and/or extra dimension) polarizations of constants 
and metric coefficients. Section V is devoted to such solutions of the 5D Einstein - Dirac equations which are 
constructed as generalizations of Taub NUT anisotropic vacuum metrics to configurations with Dirac spinor energy- 
momentum source. In Section VI we prove that the Dirac spinor field in such anisotropic spacetimes can also induce 
3D dimensional solitons which can treated as an anisotropic, in general, non-trivially typological, soliton-Dirac wave 
packet configuration propagating self-consistently in anisotropically deformed Taub NUT spacetime. Finally, in 
Section VII we conclude the work. 



II. EINSTEIN-DIRAC EQUATIONS WITH ANHOLONOMIC VARIABLES 

In this Section wc introduce an ansatz for pseudo Riemannian off-diagonal metrics and consider the anholonomic 
transforms diagonalizing such metrics. The system of field Einstein equations with the spinor matter energy- 
momentum tensor and of Dirac equations are formulated on 5D pseudo-Riemannian spacetimes constructed as a 
trivial extension by the time variable of a 4D Riemannian space (an anisotropic deformation of the Taub NUT 
instanton [^). 
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A. Ansatz for metrics 



We consider a 5D pseudo-Riemannian spacetime of signature (+,—,—,—, — ), with local coordinates 

u" = ix\y'') = (x" = = r,x^ ^ 0,y^ ^ s,y^ = p), 

- or more compactly u — {x,y) - where the Greek indices are conventionally split into two subsets a;* and t/° labeled 
respectively by Latin indices of type «, j, fc, ... ~ 0, 1, 2 and a, b, ... ~ 3, 4. The 5D (pseduo) Riemannian metric 

ds^ = gafjdu°'du^ (1) 

is given by a metric ansatz parametrized in the form 



9al3 



1 

wi /is + ni 712/14 g2 + hg + hi 102/13 7i2/i4 

W1/13 W2/13 /13 

ni/i4 n2/i4 /14 



where the coefficients are some functions of type 

9i,2 = ffi,2(a;\ a;^), ^3,4 = ^3,4(a;\ a;^, s), 
wi,2 = ^1,2(2;"^, a;^, s), 711^2 = ni^2ix^,x'^, s). 
Both the inverse matrix (metric) as well the metric (||) is off-diagonal with respect to the coordinate basis 



da 



d 
du°' 



dy"'' 



and, its dual basis, 



= du" = {d' = dx\d'' = dy"). 
The metric (|f|) with coefficients (|^) can be equivalently rewritten in the diagonal form 

5s^ = df + .gi (x) {dx^f + 52 (2:) (da;^)^ 
+/13 ((5y=^)2 + /i4 (x, s) (5^)2, 
if instead the coordinate bases and (H) we introduce the anholonomic frames (anisotropic bases) 

d 



and 



6°" EE 6u°' 

where the TV-coefficients are parametrized 

^0 = 



{5' = dx\5"' = dy" + N^{u) dx'') 



0, = wia and TVi 2 = nia 



(2) 

(3) 

(4) 
(5) 

(6) 
(7) 
(8) 



and define the associated nonlinear connection (N-connection) structure, see details in Refs |^,|l|,||,Q. The anisotropic 
frames (^ and (^) are anholonomic because, in general, they satisfy some anholonomy relations, 



with nontrivial anholonomy coefficients 



wt, = 0, w^^ = 0, = 0, w^, - 0, w:, = 0, 



(9) 



where 



is the N-connection curvature. Conventionally, the N-coefhcients decompose the spacetime variables (tensors, spinors 
and connections) into sets of mixed holonomic-anholonomic variables (coordinates) provided respectively with 'holo- 
nomic' indices of type i, j, k, ... and with 'anholonomic' indices of type a, b, c, .... Tensors, metrics and linear connections 
with coefficients defined with respect to anisotropic frames (|^) and (||) are distinguished (d) by N-coefficients into 
holonomic and anholonomic subsets and called, in brief, d-tensors, d-metrics and d-connections. 
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B. Einstein equations with anholonomic variables 
The metric (|l|) with coefficients (||) (equivalently, the d~metric (||)) is assumed to solve the 5D Einstein equations 

Ra/l - -j^dapR = uTap, (10) 

where k and T^^ are respectively the coupling constant and the energy-momentum tensor. 

The nontrivial components of the Ricci tensor for the metric (|l]) with coefficients (||) (equivalently, the d-metric 
(i)) are 

Rl^ R2= - 7; [92 



71.92 _ (ffi) 
252 2gi 



-ffi - — J> (11) 



E. /? "2 

^4i = -|^[nr+7«I], (14) 

-^42=-^ K* +7'^2] > 

where, for simphcity, the partial derivatives are denoted h* = dh/dx^, f — df /dx^ and /* = df/ds. 
The scalar curvature is computed 

R^2{Rl+ Rl) . 

In result of the obtained equalities for some Ricci and Einstein tensor components, we conclude that for the metric 
ansatz (H) the Einstein equations with matter sources are compatible if the coefficients of the energy-momentum 
d-tensor give with respect to anholonomic bases (0) and (||) satisfy the conditions 

T° = T} + Ti^,T} = = Ti,T3 = T| = T3, (15) 

and could be written in the form 

(16) 
(17) 
(18) 
(19) 



Rl- 


-KT3 


Rl- 


— kTi 















where i — 1,2 and the left parts are given by the components of the Ricci tensor (|ll[)-(|14[). 

The Einstein equations (|l0|), equivalently ([l^)-(|l9|), reduce to this system of second order partial derivation equa- 
tions: 

.. 9i92 {92? 

92 



2gi 252 

" 9l92 {9lY o 

9i - -7, 7, — = -2.9132X3, 20 

252 251 

^w, + a, = -2/i4KT3„ (22) 

2/l3 

n** + in* = — — kT4j, (23) 
/14 
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where 



hi fhl 




2 Us 




hi /h'^ 




2 Us 




(hir 


hlhl 


2hi 


2/13 ' 



"i = /^r-^(^ + ^), (24) 

a2 = - ^ ( 7^ + ? ) , (25) 



2 /l4 /l3 

and the partial derivatives are denoted, for instance, 



/3 = /ir-^-^, (26) 

lA¥-¥\ (27) 



9*2 = dg2/dx^ = dg2/dr,g[ = dgi/dx^ = dgi/O, 
hi = dhs/ds = dhs/dip (or dh^/dy^, for s = y'^). 



C. Dirac equations in anisotropic spacetimes 

The problem of definition of spinors in locally anisotropic spaces and in spaces with higher order anisotropy was 
solved in Refs. Q. In this paper we consider locally anisotropic Dirac spinors given with respect to anholonomic 
frames with associated N-connection structure on a 5D (pseudo) Riemannian space 1/(1'^'^) constructed by a direct 
time extension of a 4D Riemannian space with two holonomic and two anholonomic variables. 

Having an anisotropic d-metric 

gafiiu) = {gtj{u),hab{u)) = (l,5-(M),ft.a(u)), 

l,2;i = 0, l,2;a = 3,4, 

defined with respect to an anholonomic basis we can easily define the funfbcin (pentad) fields 

4 = f^S^ - {/i = fa = /|54, (28) 



satisfying the conditions 



9ij = flffdil and hab = faf^hab, 



gij = diag[l, —1 — 1] and hab = diag[—l, — 1]. 
For a diagonal d-metric of type (||) we have 

fi = \/\9i\^i and = \/\h^\S^, 



where 5~ and are Kronecker's symbols. 

We can also introduce the corresponding funfbein which are related with the off-diagonal metric ansatz (|^) for ga/3, 



satisfying the conditions 



and = e^a^ (29) 



7q/3 = efe^5a£ for gap = diag[l, -1, -1, -1, -1], 



ef = and ef = S^. 
The Dirac spinor fields on locally anisotropic deformations of Taub NUT spaces, 
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where / = 0, 1, are defined with respect to the 4D Euclidean tangent subspace belonging the tangent space to V^(i'2'2)^ 
The 4x4 dimensional gamma matrices 7— — [7-!- , 7- , 7- , 7- ] are defined as to satisfy the relation 



(30) 



where |7— 7- | is a symmetric commutator, g—— — (— 1, — 1, — 1, — 1), which generates a Clifford algebra distin- 
guished on two holonomic and two anholonomic directions (hereafter the primed indices will run values on the 
Euclidean and/or Riemannian, 4D component of the 5D pseudo-Riemannian spacetime). In order to extend the (pfl) 
relations for unprimed indices a, (3... we conventionally complete the set of primed gamma matrices with a matrix 7^, 
i. .e. write 7— = [7-, 7-!-, 7-, 7-, 7-] when 



{7^, 7^} = 2g^- 



The coefficients of gamma matrices can be computed with respect to coordinate bases (Q) or with respect to 
anholonomic bases (0) by using respectively the funfbein coefficients ( p8|) and (|29|), 

7"(^) = 4(^)7^ and 7^(u) = /|(z.)7^, 

were by 7"(m) we denote the curved spacetime gamma matrices and by ^^{u) we denote the gamma matrices adapted 
to the N-connection structure.. 

The covariant derivation of Dirac spinor field 5* (u) , Vq,^', can be defined with respect to a pentad decomposition 
of the off-diagonal metric (0) 



where the coefficients 



da + ^C'a^y [u) (u) 7^7^ 



(31) 



are called the rotation Ricci coefficients; the covariant derivative D-y is defined by the usual Christoffel symbols for 
the off-diagonal metric. 

We can also define an equivalent covariant derivation of the Dirac spinor field, by using pentad decompositions 

of the diagonalized d-metric (^, 



5c + \c%Au) /^(«)7^7^ 



4 2^2. 

where there are introduced N-elongated partial derivatives and the coefficients 



(32) 



are transformed into rotation Ricci d-coefhcients which together with the d-covariant derivative ' are defined by 
anholonomic pentads and anholonomic transforms of the Christoffel symbols. 

For diagonal d-metrics the funfbein coefficients can be taken in their turn in diagonal form and the corresponding 
gamma matrix 7" (u) for anisotropic curved spaces are proportional to the usual gamma matrix in flat spaces 7—. 
The Dirac equations for locally anisotropic spacetimes are written in the simplest form with respect to anholonomic 
frames, 

(i7" {u) V: - m)* = 0, (33) 
where /i is the mass constant of the Dirac field. The Dirac equations are the Euler equations for the Lagrangian 

£(1/2) (u) ^ {u) T {u) (u) (34) 

~(V;^^'+ iu))T (") * (u)] - (1-^+ (u) * (m)}, 
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where by (u) we denote the complex conjugation and transposition of the column ^ (u) 
Varying C^^^"^"^ on d-metric (34) we obtain the symmetric energy-momentum d-tensor 



(35) 



We choose such spinor field configurations in curved spacetimc as to be satisfied the conditions (|15|). 

One can introduce similar formulas to (^3|)-(|3^) for spacetimes provided with off-diagonal metrics with respect to 
holonomic frames by changing of operators 7^ (u) 7^ (u) and V/3 —^ V/3. 



III. ANISOTROPIC TAUB NUT - DIRAC SPINOR SOLUTIONS 

By straightforward calculations we can verify that because the conditions D^^ fa — ^ ^^''^ satisfied the Ricci rotation 



coefficients vanishes, 



4*^7 (^) - and V:^^' - Sa'f, 



and the anisotropic Dirac equations ( |33| ) transform into 

(^7" (u) Sa - Ai)* = 0. (36) 

Further simplifications are possible for Dirac fields depending only on coordinates {t, = r, = 0), i. e. ^ — ^'(a;*'') 
when the equation (^6|) transforms into 

+ Z7i— + ij^^=d2 - m)* - 0. 
Vl5i| ^152 

The equation (|3^) simplifies substantially in ^-coordinates 
defined as to be satisfied the conditions 

= ^r=f^i ^'^'^ 4p2 = ^r=f^2 (37) 
9^ \/\m\ J\g^\ 



We get 

(-*7o| + + ^l2-^ - A*)*(i, C\ C') = 0. (38) 

The equation ( p8[ ) describes the wave function of a Dirac particle of mass fi propagating in a three dimensional 
Minkowski flat plane which is imbedded as an anisotropic distribution into a 5D pseudo-Riemannian spacetime. 
The solution * = *(t, C\ C^) of (|3) can be written 



*(+)(C) = exp[-i(fcoi + fciC' +fc2C')]'^°(fc) 
for positive energy; 

v&(-)(C) - exp[z(M + fciC'+fc2C')]x"(fc) 
for negative energy, 



with the condition that fco is identified with the positive energy and ^'^{k) and 'X^(k) are constant bispinors. To satisfy 
the Klein-Gordon equation we must have 

Jxi — ^2 — ' 
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The Dirac equations implies 



{a'h - ^i)^''{k) and [a'h + Ai)x°(fc), 



where a^{i — 0, 1, 2) are Pauh matrices corresponding to a reahzation of gamma matrices as to a form of sphtting to 
usual Pauli equations for the bispinors and 

In the rest frame for the horizontal plane parametrized by coordinates ( — {i, C^iC^} there are four independent 
solutions of the Dirac equations, 



<y9^2)(M,0) 



In order to satisfy the conditions ( p^ ) for compatibility of the equations (p0|)-(p3[) we must consider wave packets of 
type (for simplicity, we can use only superpositions of positive energy solutions) 



*(+)(() = 



X fc(p,W)^'"l(fc)exp[-ife.C] (39) 

[a] = l,2,3 



when the coefficients b{p, [a]) define a current (the group velocity) 



with Ip^I ^ and the energy-momentum d~tensor ( |35| ) has the next nonrivial coefficients 

Tj] = 2T(C\C') =A:o*+7o*, 

T} = -ki^+-fi^, Tl = -/c2*+72* (40) 

where the holonomic coordinates can be reexpressed Q — C{x^)- We must take two or more waves in the packet and 
choose such coefficients h{p, [a]), satisfying corresponding algebraic equations, as to have in ( ^o|) the equalities 

T\ = Tl = T{C\e) = r{x\x^), (41) 

required by the conditions (^5|). 



IV. TAUB NUT SOLUTIONS WITH GENERIC LOCAL ANISOTROPY 

The Kaluza-Klein monopole was obtained by embedding the Taub-NUT gravitational instanton into five- 
dimensional theory, adding the time coordinate in a trivial way. There are anisotropic variants of such solutions j^] 
when anisotropies are modelled by effective polarizations of the induced magnetic field. The aim of this Section is to 
analyze such Taub-NUT solutions for both cases of locally isotropic and locally anisotropic configurations. 



A. A conformal transform of the Taub NUT metric 



We consider the Taub NUT solutions and introduce a conformal transformation and a such redefinition of variables 
which will be useful for further generalizations to anisotropic vacuum solutions. 
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1. The Taub NUT solution 



This locally isotropic solution of the 5D vacuum Einstein equations is expressed by the line element 

ds%D) = dt" + dsl^jjy (42) 
ds^4^) = -V-^{dr^ + r'^de'^ + sin^ Odip^) 
-V{dx^ + A.dx'f 

where 

V = 1 H ,?7io = const. 

r 

The functions Ai are static ones associated to the electromagnetic potential, 

Ar = 0, = 0, A^ = 4mo (1 - cos 9) 

resulting into "pure" magnetic field 

B = rot A = . (43) 

of a Euclidean instanton; 1^ is the spherical coordinate's unity vector. The spacetime defined by (^2|) has the global 
symmetry of the group Gg = 5*0(3) ® Ui{l) ® Tt(l) since the line element is invariant under the global rotations of 
the Cartesian space coordinates and and t translations of the Abehan groups C/4(l) and Tt{l) respectively. We 
note that the Ua{\^) symmetry eliminates the so called NUT singularity if y'^ has the period Anmo- 



2. Conformally transformed Taub NUT metrics 



With the aim to construct anisotropic generalizations it is more convenient to introduce a new 5th coordinate, 

y^^^ = y^- j fi-\d,^)dae,^), (44) 

with the property that 



which holds for 



when 



and, for instance. 



d<^ + 4too(1 — cos6)d9 = dy'^ + 4r7io(l — C089)dip, 



d^ = fiie, ip)di, ~ y^) - l^rffl + ^d^, 



— = 4mo(l - cos 9)fi, 

— = -4mo(l - cos0)/Lt, 



^ = (1 — cos 9) ^ exp[9 — ip] 



The changing of coordinate ( [44[ ) describe a reorientation of the 5th coordinate in a such way as we could have only 
one nonvanishing component of the electromagnetic potential 

Ag = 4mo (1 - cos 6*) . 

The next step is to perform a conformal transform, 
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and to consider the 5D metric 

ds%D) = dt^ + dsf^D)- (45) 
= ~{dr^ + r^de^) - sin^ Odip^ 

~V^{dC + Aed9f, 

(not being an exact solution of the Einstein equations) which wiU transform into some exact solutions after corre- 
sponding anholonomic transforms. 

Here, we emphasize that we chose the variant of transformation of a locally isotropic non-Einsteinian metrics into an 
anisotropic one solving the vacuum Einstein equations in order to illustrate a more simple procedure of construction 
of 5D vacuum metrics with generic local anisotropy. As a metter of principle we could remove vacuum isotropic 
solutions into vacuum anisotropic ones, but the formula in this case would became very combersome.. The fact of 
selection as an isotropic 4D Riemannian background just the metric from the linear interval ^5^4^^ can be treated as 
a conformal transformation of an instanton solution which is anisotropically deformed and put trivially (by extension 
to the time like coordinate) into a 5D metric as to generate a locally isotropic vacuum gravitational field. 



B. Anisotropic Taub NUT solutions with magnetic polarization 

We outline two classes of exact solutions of 5D vacuum Einstein equations with generic anisotropics (see details in 
Ref. [^) which will be extended to configurations with spinor matter field source. 



1. Solutions with angular polarization 

The ansatz for a d-metric (^), with a distinguished anisotropic dependence on the angular coordinate ip, when 
s = if, is taken in the form 

Ss^ = dt^ - SsJ 



Ssljj) = -{dr^ + r^d9^) - sin^ 6*^ 



'(4D) 

'■(4D) - '"""^ 

(5<; = d"; + n2{9, (p)d9, 



where the values f^K^, (p) (we use non-negative values 77I not changing the signature of metrics) and n2{9, (p) must be 
found as to satisfy the vacuum Einstein equations in the form (|20|)-(|23|). We can verify that the data 



x° ^t,x^ = r, = 9,y^ ^s^ ip, = (46) 
go = 1,51 = -1,32 = -r^-.hs = -r^ sin^ 6*, 
hi = (r) 7?2^) , 772^) = [1 + zu{r, 9)^f,w, - 0; 

no,i = 0; 712 = n2[o] {r, 9) + n2[i] (r, 9) /[I + w{r, 9)Lpf . 



give an exact solution. If we impose the condition to obtain in the locally isotropic limit just the metric (45), we have 
to choose the arbitrary functions from the general solution of (|2^) as to have 

77^^) = [1 + uj{r, 9)^f ^ 1 for Tu{r, 9)^ ^ 0. 
For simplicity, we can analyze only angular anisotropics with vj = voij))^ when 

riU=ril^){(^.^) = [l + ^{e)p>Y. 

In the locally isotropic limit of the solution for ?i2 (f, 9, ip), when vjip — > 0, we could obtain the particular magnetic 
configuration contained in the metric ( p5| ) if we impose the condition that 

«2[o] {r, 9) + 772(1] (f , 9) = Ag = 4777-0 (1 - cos 9) , 
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which defines only one function from two unknown values n2[o] {t, 0) and n2[i] (r, 9) . This could have a corresponding 
physical motivation. From the usual Kaluza-Klein procedure we induce the 4D gravitational field (metric) and 4D 
electromagnetic field (potentials Ai), which satisfy the Maxwell equations in 4D pseudo-Riemannian spacetime. For 
the case of spherical, locally isotropic, symmetries the Maxwell equations can be written for vacuum magnetic fields 
without any polarizations. When we introduce into consideration anholonomic constraints and locally anisotropic 
gravitational configurations the effective magnetic field could be effectively renormalized by higher dimension grav- 
itational field. This effect, for some classes of anisotropics, can be modeled by considering that the constant mo is 
polarized, 

mo -> m (r, 9, ip) = mo?]™ (r, 9, ip) 

for the electromagnetic potential and resulting magnetic field. For "pure" angular anisotropics we write that 

712 {9, ip) = n2[o] {0) + ri2[i] {0) /[I + ^{9)^? 
= 4mo?7„ {9,Lp) (1 - cos 6*) , 

for 

77^^) (0, p) = 772^)[o] {9) + 772^)[i] {9) /[I + uj{9)p]\ 
This could result in a constant angular renormalization even zu{9)ip —>■ 0. 



2. Solutions with extra-dimension induced polarization 

Another class of solutions is constructed if we consider a d-metric of the type (^), when s = i;, with anisotropic 
dependence on the 5th coordinate <j, 

Ssf^jj) = -{dr^ + r^d9^) ~ sin^ 9dLp^ 
-V'ir)ril^i9,,)Se, 
(5? = d^ + W3{9,<;)d9, 

where, for simplicity, we omit possible anisotropics on variable r, i. e. we state that ?7(^-) and W2 are not functions of 
r. 

The data for a such solution are 

x° =t,x^= r, = 9,y^ = s = ^, y* = ^, (47) 
go = l,gi = -1,52 = -r'^.hi = -r^ sin^ 6*, 
7i3 = "1/2 {r)ijf^),'qf^) = 77^^) (r,0,<;), 710,1 = 0; 
7^0,1 = 0,W2 = 4,rnori„r (0, ^) (1 - cos 6*) ,7io = 0, 

"1,2 = »T-i,2[o] {r, 0) + n-i,2[i] (»', 0) j ^)*, 

where the function Ty^,;) — 77(<;) (r, 9, <;) is an arbitrary one as follow for the case /i| — 0, for angular polarizations we 
state, for simplicity, that 77(5) does not depend on r, i. e. Ty^,;) — 77(^)(0, <j). We chose the coefficient 

W4 = 4mo7ym {9, ^) (1 - cos 6*) 

as to have compatibility with the locally isotropic limit when W2 — Ag with a "polarization" effect modeled by 
r]m {9, <;) , which could have a constant component rjm — Vmlo] = const for small anisotropics. In the simplest cases we 
can fix the conditions 7ii_2[o,i] i'r,9) = 0. All functions rj'^^^rjm and t^i, 2(0,1] can be treated as some possible induced 
higher dimensional polarizations. 
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V. ANISOTROPIC TAUB NUT-DIRAC FIELDS 



In this Section we construct two new classes of solutions of the 5D Einstein-Dirac fields in a manner as to extend 
the locally anisotropic Taub NUT metrics defined by data (^) and ( ^ ) as to be solutions of the Einstein equations 
(^)^(p3|) with a nonvanishing diagonal energy momentum d-tensor 

= {2T{r,0),T(r,e),T{r,e),O,O} 

for a Dirac wave packet satisfying the conditions (^) and (p]). 



A. Dirac fields and angular polarizations 

In order to generate from the data ( p6| ) a new solution with Dirac spinor matter field we consider instead of a linear 
dependence of polarization, ryj^) ^ [l + vu (r, 9) ip], an arbitrary function rj^^-) (r, 0, ip) for which hi = V^{T^)V(ip) ¥') 
is an exact solution of the equation ( pl| ) with Ti = T (r, 9) . With respect to the variable rj^^-^ (r, 9, ip) this component 
of the Einstein equations becomes linear 



(48) 



which is a second order linear differential equation on variable (p with parametric dependencies of the coefficient 
sin^ 9T on coordinates (r, 9) . The solution of equation (Es) is to be found following the method outlined in Ref. 



Ci (r, 9) coah[ipr sin6'^|T (r, 
T(r,0) <0; 

Ci ir,9)+C2ir,9)^,T{r, 



C2{r,9)], 



0) = O; 

= Ci ir,9) cos[^rsm9./T^+C2 (r, 9)], 
T{r,9)>0, 



(49) 
(50) 

(51) 



where Ci^2 d) are some functions to be defined from some boundary conditions. The first solution ([49|), for negative 
densities of energy should be excluded as unphysical, the second solution (50) is just that from (46) for the vacuum 
case. A new interesting physical situation is described by the solution ( ^l| ) when we obtain a Taub NUT anisotropic 
metric with periodic anisotropic dependencies on the angle ip where the periodicity could variate on coordinates (r, 9) 
as it is defined by the energy density T (r, 9) . For simplicity, we can consider a package of spinor waves with constant 
value of T = Tq and fix some boundary and coordinate conditions when Ci^2 = Ci 2[o] are constant. This type of 
anisotropic Taub NUT solutions are described by a d-metric coefficient 



hi = y'^{r)C^o] cos^ [ipr sin 9 y 



'To + C2 



(52) 



Putting this value into the formulas (p^), ( |2^ ) and (|2^) for coefficients in equations ( |22| ) we can express ai^2 — 
cti,2[h3, /i4, To] and /? = Pihs, hi, To] (we omit these rather simple but cumbersome formulas) and in consequence we 
can define the values w\^2 by solving linear algebraic equations: 

wi,2 [r, 9, (f) = ai^2 {r, 9, ip) /(3 {r, 9, ip) . 

Having defined the values ( p^ ) it is a simple task of two integrations on ip in order to define 

1 + cos K 1 1 



"2 = "2[o] {1; 0) In ■ 
+"2[i] [r, 9) 



1 — cos k 1 — cos k 1 — sin k 



(53) 



were 



— (/?r sinf^yTo + C; 



2[0], 

n2[o,i] {r,9) are some arbitrary functions to be defined by boundary conditions. We put no,i = to obtain in the 
vacuum hmit the solution (B^). 
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Finally, we can summarize the data defining an exact solution for an anisotropic (on angle ip) Dirac wave packet - 
Taub NUT configuration: 





— i — r — f) ifl — a — 

— t J Jl - — - / J Jb - — - u , y — o — 


4 

2/ — 


50 


= 1,51 = -1,52 = -7-^,/j3 


-r^ sin^ 9, 


hi 




cos k{r, 6, (fi), 




= 0, no,i = 0, ^2 = n2 [r, 6, k 


(r,6',(p)) see (53), 




= ¥+^ {(^^-^x^x^)) see (I 




T 


= T{C'^{x\x^)) see®. 





This solution will be extended to additional soliton anisotropic configurations in the next Section. 



B. Dirac fields and extra dimension polarizations 

Now we consider a generalization of the data ( p7[ ) for generation of a new solution, with generic local anisotropy on 
extra dimension 5th coordinate, of the Einstein - Dirac equations. Following the equation (g2|) we conclude that there 
are not nonvacuum solutions of the Einstein equations (with T 7^ 0) if /i| = which impose the condition T = for 
^3, ^4 7^ 0- So, we have to consider that the d~metric component ft-4 = —r^ sin^ 9 from the data ( ^7| ) is generalized to 
a function hi (r, 6*, <;) satisfying a second order nonlinear differential equation on variable <: with coefficients depending 
parametrically on coordinates {r,9) . The equation (i. e. (p2|)) can be linearized (see Ref. if we introduce a new 
variable hi — h^, 

h** - ^h* + hzTh = 0, 

which, in its turn, can be transformed : 

a) to a Riccati form if we introduce a new variable w, for which h = v* /v, 



v*+v^- -^v + hsT = 0; 



b) to the so-called normal form ||T^ , 



(55) 



A** +/A = 0, 



(56) 



obtained by a redefinition of variables like 



A = /lexp 



1 fK/ 



— h h. 



^1/4 



where 



16/13 



We can construct explicit series and/or numeric solutions (for instance, by using Mathematica or Maple programs) 
of both type of equations ( p5| ) and normal ( p6| ) for some stated boundary conditions and type of polarization of the 
coefficient h^, (r, 9, <j) = V"^ (nvf^) ^, ^) ^'^d, in consequence, to construct different classes of solutions for hi (r, 9, (^) . 
In order to have compatibility with the data ( ^7| ) we must take /14 in the form 



hi ir,9,(i) 



-r^ sin^ I 



+ /i4(<;) ir,9,<;) . 



where hif^\ (r, 9, <;) vanishes for T — > 0. 



Having defined a value of hi (r, 9, <j) we can compute the coefficients (p^, ( [25| ) and ( p6| ) and find from the equations 



wi,2 (f , 9, <;) = ai,2 (r, 9,<,) /(3 (r, 9, <,) . 
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From the equations (p3|), after two integrations on variable <; one obtains the values of ni.2 {r, 0, <;) . Two integrations 
of equations (^3|) define 



where 



ni{r,9,i;) = ni[o]{r,9) dz dsP{r,9, s) + ni[i]{r,t 
Jo Jo 



1 h* h* 

2 /i3 n.4 



and the functions ni[Q]{r,9) and ni[^{r,9) on (r, 6*) have to be defined by solving the Cauchy problem. The boundary 
conditions of both type of coefhcients Wi 2 and ni 2 should be expressed in some forms transforming into corresponding 
values for the data ( p7| ) if the source T ^ 0. We omit explicit formulas for exact Einstein-Dirac solutions with 
polarizations because their forms depend very strongly on the type of polarizations and vacuum solutions. 

VI. ANHOLONOMIC SOLITON-TAUB NUT-DIRAC FIELDS 

In the next subsections we analyze two explicit examples when the spinor field induces two dimensional, depending 
on three variables, solitonic anisotropies. 

A. Kadomtsev— Petviashvili type solitons 

By straightforward verification we conclude that the d-metric component hi(r,9,s) could be a solution of 
Kadomtsev-Petviashvili (KdP) equation jisj (the first methods of integration of 2+1 dimensional soliton equations 
where developed by Dryuma and Zakharov and Shabat pO[|) 

K* + e(hi + Qh^h'^ + h'A ' = 0, e = ±1, (57) 



if the component h^{r,9,s) satisfies the Bernoulli equations |17 

hl + Y (r, 9, s) (hsf + F, (r, 9, s) = 0, (58) 

where, for /i| ^ 0, 

Y{r,9,s)^KT^^, (59) 



F, (r, 0, s) = -1 + -i + 6/i4/i4 + h'l' 



and 

The three dimensional integral variety of (pq) is defined by formulas 



where 



h-' {r, 9, s) = h-^l^ (r, 9) {x\ s) x J ^^^^ds, 
E^ {r, 9, s) — cxp J Fe (r, 9, s) ds 

and /i3(a;) (r, 9) is a nonvanishing function. 

In the vacuum case Y (r, 9,s) = and we can write the integral variety of (Fsl) 



/i^™^) (r, 0, s) = h^llf (r, 9) exp - J F, (r, 9, s) ds 



We conclude that a solution of KdP equation (^) could be generated by a non-perturbative component ft-4(r, 9, s) 
of a diagonal h-metric if the second component (r, 9, s) is a solution of Bernoulli equations (|5^) with coefficients 
determined both by and its partial derivatives and by the T\ component of the energy-momentum d-tensor (see 
(p])). The parameters (coefficients) of (2+1) dimensional KdP solitons are induced by gravity and spinor constants 
and spinor waves coefficients defining locally anisotropic interactions of packets of Dirac's spinor waves. 
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B. (2+1) sine— Gordon type solitons 



In a similar manner we can prove that solutions /i4(r, 6, s) of (2+1) sine-Gordon equation (see, for instance, pl|-p3[) 

^4* + — h4 — sin(ft,4) 
also induce solutions for /13 (r, 0, s) following from the Bernoulli equation 

h; + KEir, 9)^{h3)^ + F (r, 9, s) ^0,hl^ 0, 
"■4 

where 

^,4 — ft,4 — sin(/i4) . 

The general solutions (with energy-momentum sources and in vacuum cases) are constructed by a corresponding 
redefinition of coefficients in the formulas from the previous subsection. We note that we can consider both type of 
anisotropic solitonic polarizations, depending on angular variable ip or on extra dimension coordinate <j. Such classes 
of solutions of the Einstein-Dirac equations describe three dimensional spinor wave packets induced and moving self- 
consistently on solitonic gravitational locally anisotropic configurations. In a similar manner, we can consider Dirac 
wave packets generating and propagating on locally anisotropic black hole (with rotation ellipsoid horizons), black 
tori, anisotropic disk and two or three dimensional black hole anisotropic gravitational structures [|lj. Finally, we note 
that such gravitational solitons are induced by Dirac field matter sources and are different from those soliton solutions 
of vacuum Einstein equations originally considered by Belinski and Zakharov pj] . 



F(r,( 



n.4 



h* 



VII. CONCLUSIONS 



We have argued that the anholonomic frame method can be applied for construction on new classes of Einstein-Dirac 
equations in five dimensional (5D) spacetimes. Subject to a form of metric ansatz with dependencies of coefficients on 
two holonomic and one anholonomic variables we obtained a very simplified form of field equations which admit exact 
solutions. We have identified two classes of solutions describing Taub NUT like metrics with anisotropic dependencies 
on angular parameter or on the fifth coordinate. We have shown that both classes of anisotropic vacuum solutions 
can be generalized to matter sources with the energy-momentum tensor defined by some wave packets of Dirac fields. 
Although the Dirac equation is a quantum one, in the quasi-classical approximation we can consider such spinor fields 
as some spinor waves propagating in a three dimensional Minkowski plane which is imbedded in a self-consistent 
manner in a Taub-NUT anisotropic spacetime. At the classical level it should be emphasized that the results of this 
paper are very general in nature, depending in a crucial way only on the locally Lorentzian nature of 5D spacetime and 
on the supposition that this spacetime is constructed as a trivial time extension of 4D spacetimes. We have proved that 
the new classes of solutions admit generalizations to nontrivial topological configurations of 3D dimensional solitons 
(induced by anisotropic spinor matter) defined as solutions Kadomtsev-Petviashvili or sine-Gordon equations. 
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